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I. INTRODUCTION 



In previous papers we analyze a class of models of D-dimensional spherically 

symmetric random walks, where D is not restricted to integer values. In this paper we 
extend these models to allow for the creation and annihilation of random walkers. We 
demonstrate that these extended models exhibit critical behavior as a function of the birth 
rate of walkers. The critical coefficients depend on the value of the dimension D. In the 
next paper in this series we apply these ideas to the study polymer growth in D dimensions 
in the vicinity of a hyperspherical adsorbing boundary. 

The random walks in Refs. take place on an infinite set of regions labeled by the 

integer n, n = 1, 2, 3, .... If the random walker is in region n at time t, then at time 
t + 1 the walker must move outward to region n + 1 with probability Pout{n) or inward to 
region n — 1 with probability Pi^{n), where 

PoutW+PinH = l (1.1) 

SO that probability is conserved. [We take Pout(l) = 1 and Pin(l) = to enforce the 
requirement that a walker in the central region (n = 1) must move outward at the next 
step.] Let Cn,t;m represent the probability that a random walker who begins in the mth 
region at t = will be in the nth region at time t. The probability Cn,t;m then satisfies the 
difference equation 

r Pin(ra + l)C„+i,t_i;„ + Pout('^ - l)C„_i,t-i;m {n > 2), 

Cn,t;m = \ (1-2) 

[ Pi„(2)C2,t-i;™ {n = 1) 

and the initial condition 

Cn,0;m — ^n,m- (-'-•3) 

To formulate a model of spherically symmetric random walks in D-dimensional space we 
take region n to be the volume bounded by two concentric D-dimensional hyperspherical 
surfaces of radii Rn-i and P„. In Ref. [|l| we take the probabilities of moving out or in to be 



in proportion to the hyperspherical surface areas bounding region n. Let Sd{R) represent 
the surface area of a D-dimensional hypersphere 



Then, for n > 1, 



P (ti) - Spj^n) _ Rn ^ 



and 



So{Rn)+So{Rn^i) R^-' + R^ll' 
As we discussed earher, for the special case n = 1 we define 

Pout(l) = l and Pin(l)=0. (1.6) 



The choices in Eqs. (p..4|)-( pr6D satisfy the requirement that probabihty be conserved because 
they obey Eq. ([TTI ). 

For dimensions other than D = 1 and D = 2, when we substitute Eqs. ( p..4D -( pr^ into 
Eq. ( |1.2D and take Rn = n, we obtain a difference equation that cannot be solved in closed 
form. Thus, in Ref. we proposed that the probabilities in Eqs. ( |1.4|) - (|1.6|) be replaced by 
bilinear functions of n, which are are uniformly good approximations to Pout(^) and Pin{n) 
in the range D > when Rn = n: 

PouM) = !!\^~\ and Pin(n)= ^"^ (1.7) 
2n + D — 6 2n + D — 6 

Now, the difference equation initial- value problem ( |1.2|) and ( p..3|) for the probabilities Cn,t;m 
can be solved in closed form: 

Cn,t;m = 2D-ir2(£) /2)T{m + D -2) 1-1 ^ ^^''^^ "^^^^^^^i-i \x)cL-i (1.8) 

where (x) is a Gegenbauer polynomial ^ . From the solution in Eq. ( |1.8D one can obtain 
closed- form expressions for spatial and temporal moments of the random walk Q. 



In this paper we generalize the difference equation ( p..2|) to include the possibility of 
creation and annihilation of random walkers. We allow random walkers to give birth in 
region 1 with birth rate a and to die in all other regions with uniform death rate z. Birth 
rates and death rates are properties of populations rather than of single individuals. Thus, 
rather than solving Eq. ( |1.2[ ) as an initial-value problem for a single random walker who 
starts in region m, we are going to study a large population of random walkers, all of 
whom obey this difference equation. We represent this population of random walkers by a 
distribution Gn,t, which denotes the number of random walkers in region n at time t. The 
distribution Gn^t satisfies the same recursion relation as Cn,t;m except for the factors of a 
and z: 

' zPin{n + l)Gn+i,t-i + zPout{n - l)G'n-i,t-i (n > 3), 

Gn,t = I zPi^i3)G3,t-i + aGi,t-i {n = 2), (1.9) 

^ zP,^{2)G2,t-i (n = 1), 

where we have set Pout(l) = 1- Note that the function Gn,t must be positive for all n 

and t. Aside from the requirement that Gnfl, the initial distribution of random walkers, be 

normalizable it is arbitrary. We are not concerned with the detailed structure of the initial 

distribution; rather, we are interested in the asymptotic behavior of distributions as t — cxd. 

The specific choice of the initial distribution is unimportant because a random walk is a 

diffusive (dissipative) process and details of Gnfi are irretrievably lost as time evolves; all 

initial distributions lead to the same large-time behavior. This behavior is determined by 

the details of the random walk process itself. 

In this model random walkers are created or destroyed at a given site in proportion 

to the number of walkers at that site, where a and z are the constants of proportionality. 

Technically speaking, a acts as a birth rate if a > 1; if a < 1, it is really a death rate. A 

similar interpretation applies to z. We are particularly interested in steady-state solutions 

of Eq. ( |1.9| ); the existence of such solutions imposes a relationship between the birth rate 

and the death rate. 

In Sec. |l| we perform numerical and analytical studies of Eq. ( |1.9|) for arbitrary choice 



of -Pout('^) and Pout('^)- We study two quantities: Nt, the total number of random walkers 
at time t, and Ft, the fraction of random walkers in region 1 at time t. We show that the 
(physically relevant) positive quadrant of the (a, z) plane is partitioned by two intersecting 
critical curves into four regions (or phases) characterized by the behavior of Nt and Ft 
as t ^ oo. (Some detailed asymptotic studies of the large-t behavior are given in the 
Appendix.) The intersection of these two critical curves is the critical point (oc-Zc) for 
steady-state distributions of random walkers. Although the location of this critical point is 
a function of the dimension D, the qualitative features of this phase diagram are generic. 

We obtain a solvable model of random walks for arbitrary dimension D by using the 
uniform approximations for Pout(^) and Pin{n) proposed in Ref. and given in Eq. { \L.7\) . 
In Sec. m we analyze this model for the simple case D = 1. We analytically determine 
the features of the (a, z) phase diagram and verify that random walk distributions exhibit 
critical behavior. Then, in Sec. |^ we use generating-function methods to solve Eq. ( |1.9|) 
for arbitrary D. We cannot perform a global analysis of the (a, z) phase diagram as in the 
case D = 1 but we can perform a local analysis in the vicinity of the critical point. From 
this local analysis we show that a second-order phase transition occurs at the critical point. 

Specifically, we show that near the critical birth rate Cc, the steady state distribution 
fraction F{a) = limt^oo-^t behaves like 

F{a) C{D){a-a^y (a ^ a+, D 7^ 2, 4), (1.10) 

where the multiplicative constant C{D) depends on the dimension D. The critical exponent 
u also depends on the dimension D: 

^ (0<D<2), 



z/ = < 



2-D 
2 



(2<D<4), (1.11) 



D-2 

I 1 {D> 4). 

There is no critical exponent for the special cases D = 2,4; instead, we find that as a — 
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In general, formulating simplified D-dimensional statistical models is useful for under- 
standing aspects of critical phenomena exhibited actual physical systems. Indeed, any solv- 
able statistical model that exhibits nontrivial critical behavior is worthy of study 0. In the 
next paper in this series we apply the results of this paper to the study of polymer growth 
in D dimensions. There we extend to arbitrary dimension the earlier results for D = 1 
and D = 2 

II. RANDOM WALKS WITH BIRTH AND DEATH 

In this section we discuss the general properties of the spherically symmetric random 
walk model defined by Eq. (|1.9|) in which random walkers may be created and annihilated. 
Let 

oo 

iVt = E Gn,t (2.1) 

n=l 

be the total number of random walkers at time t. We restrict our attention to initial 
distributions for which Nq is finite so that A^^^ is finite for all t. Let 

Ft = GijNt (2.2) 

represent the fraction of all random walkers in region 1 at time t. 

Numerical ||^ and analytical studies of the quantities Nf and as t — > oo reveal that, 
independent of the initial distribution of random walkers, the asymptotic behaviors of A''^ 
and Ft are determined by the values of a and z. Specifically, we obtain the generic result 
that for any value of D the positive quadrant of the (a, z) plane is partitioned into four 
distinct regions by two boundary curves as shown in Fig. 

One of the boundary curves, which we have labeled Bi in Fig. is a straight line passing 
through the origin. To the left of Bi we find that Ft vanishes as t ^ oo; to the right of Bi 
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we find that Ft approaches a positive finite value as t — > oo. For D < 2 the equation for 
the boundary fine Bi is z = a; as D increases beyond 2 the boundary fine remains straight 
but the slope of Bi begins to decrease with increasing D. As we will see, the transition that 
occurs at D = 2 is a refiection of Polya's theorem [|1^], which states that when D > 2 the 
probability of an individual random walker visiting region 1 more than once is less than 
unity. 

The second boundary curve shown in Fig. |l| is labeled B2. This curve consists of two 
parts: The first part is a straight-line segment, z = 1, extending from the z axis to the 
boundary line Bi. This line segment connects to the second part, which is a curve that 
approaches z = as a oo. The equation describing the second part depends on D. [For 
D = 1 this curve is given hj z = 2a/ [o? + 1) (a > 1), as we will show in Sec. |T|.] Above 
the boundary B2 we find that Nt 00 as t ^ 00; below B2 we find that A^^ — as t — 00. 
On B2 the total number of walkers approaches a finite distribution N{a) as t ^ 00. On the 
curved portion of B2 the function N{a) is positive; on the straight-line portion of B2 the 
function N{a) is positive for D > 2 while N{a) = for D < 2. This transition at -D = 2 is 
yet another manifestation of Polya's theorem. 

While detailed studies of the large-t asymptotic behavior of Nt and Ft are given in 
the Appendix, many of the qualitative features of Fig. [l| can be derived directly from an 
analysis of Eq. (|1.9| ). To determine the boundary line Bi we introduce a change of variable 
in Eq. (|r|): 



Gn,t — Z Hn^f (2-3) 

The distribution if„ ^ satisfies the recursion relation for a D-dimensional spherically sym- 
metric random walk with a birth rate a/z m region 1 and no births or deaths occurring in 
any other region: 

' Pin(n + l)Hn+l,t-l + Pont{n - l)Hn^l,t~l > 3), 



n,t 



Pin{^)Hs,t-i + fi/i,i-i {n = 2), (2.4) 

[PU2)H2,t-i in = l). 



Let Ili{D) denote the probability that a random walker in region 1 will eventually return 
to region 1. Suppose a random walk satisfying Eq. ( p.4| ) begins at t = 0. Of the Hi q walkers 
who begin in region 1, only a fraction Ili{D) of them will eventually return to region 1 
to give birth to new walkers at the rate a/z. Of these new walkers, again only a fraction 
Ili{D) of them will return to region 1 to give birth again, and so on. Hence, to find the 
total number of random walkers who are ever born we must sum a geometric series whose 
geometric ratio is the quantity aIli{D)/z. If this quantity is less than 1, the geometric series 
converges and the total number of random walkers ever born is finite. As time t increases, 
the random walkers diffuse away from region 1. Thus, the ratio 

* Y^OO Y^OO TT 

vanishes as t — oo. On the other hand, if the quantity aIli{D)/z is greater than 1, both 
Hi t and I]^o Hn,t diverge at the same rate and the ratio Ft approaches a nonzero limit 
(that lies between and 1) as t ^ oo. 

The transition between Ft ^ and Ft — > finite limit occurs on the line 

z = ani{D). (2.5) 

This is the equation of the boundary line Bi. Polya's theorem states that for any random 
walk Ili{D) = 1 when D < 2 and Ili{D) < 1 when D > 2. This theorem explains the 
transition in the slope of the line Bi at D = 2. In the spherically symmetric random walk 
model discussed in Ref. where Pout('^) and Pin{n) are given in Eqs. ( |1.4| - |T^ ), it was 
shown that 

Ui{D) = l-l/an-l) {D>2) (2.6) 

(C is the Riemann Zeta function); in the random walk model discussed in Ref. 0, where 
-Pout(^) and Pin{n) are given in Eq. (|1.7D, it was shown that 

Ui(D) = l/{D-l) {D>2). (2.7) 
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Numerical computation confirms the slope of the boundary line Bi for both models (see 

Figs, m- 

The shape of the curved part of the boundary B2 in Fig. ^ depends on the dimension 
D and on the choice of the functions Pout('^) and Pin(n). It is not universal. However, the 
straight-line portion of the boundary B2 is universal and is easy to understand for any D. 
Points (a, z) such a < ac and z near 1 lie to the left of Bi. Thus, Ft, the fraction of random 
walkers in region 1, becomes vanishingly small as t — 00. Hence, the effect of the birth rate 
a on the total number of walkers is negligible. The growth or decay of the total number of 
walkers only depends on the magnitude of z; if z < 1 then Nt —>■ a.s t —>■ 00, and if z > 1 
the — > 00 as t — s> CX3. 

On the straight-line portion of the curve B2, where a < a^. and z = 1, the limiting value 
of Nf depends on the dimension D. If D < 2 then = 1. Thus, on this portion of B2 a 
fraction 1 — a of random walkers who arrive in region 1 at a given time step must die at 
the next time step. But by Polya's theorem all random walkers visit region 1 repeatedly. 
Hence, the total number of random walkers Nt must vanish as t ^ 00. On the other hand, 
if D > 2 we have Ili{D) < 1. Thus, the fraction 1 — Hi(D) of random walkers who originate 
in region 1 never return to region 1. Thus, these random walkers never die because z = 1. 
Hence, Nt approaches a finite positive number as t 00. 

We find numerically that as we cross the boundary line Bi, the limiting value of the 
function as t — 00 is continuous. We are particularly interested in crossing from one 
side of Bi to the other along the boundary curve B2 that divides the upper region, where 
Nt — 00, and the lower region, where A^^ — > as t — 00. We focus on this curve B2 because 
it is only on this curve that a steady state is reached as t 00. Along this boundary 
curve the limiting value of Ft undergoes a second-order phase transition at the critical point 
(oc, Zc), which is situated at the intersection of Bi and i?2. On the curve B2 when a < ac the 
limiting value of Ft is (even though the limiting value of Nt may be 0), and when a > ac 
the limiting values of both A^^^ and Ft on the boundary curve B2 are finite positive numbers. 
The curved portion of B2 is in fact the locus of all points in the positive quadrant of the 
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(a, z) plane for which the hmiting values of both Nt and Ft as t oo are finite and nonzero. 

The interpretation of lim^^oo being finite and nonzero is that the distribution ^ 
approaches a steady-state. In such a steady-state there is a balance between random walkers 
being created in region 1 and annihilated in all other regions. This steady-state solution 
can be obtained by solving a discrete eigenvalue problem. 

Steady-state distributions are special cases of shape-independent distributions; that is, 
distributions that do not change shape as they evolve in time. For such distributions 
Gn,t/Gm,t is independent of t for all n and m so that the relative number of walkers in 
region n is a time-independent fraction of the total number of walkers. The time depen- 
dence of such distributions is very simple: 

Gn,t = 9nX'- (2.8) 

The distribution satisfies the discrete eigenvalue problem 

' Pin(n + l)zgn+i + Pontin - l)zgn-i {n > 3), 

Pin{3)zg3 + agi iri = 2), (2.9) 

[PU2)zg2 (n = l), 

which is obtained by substituting Gn,t in Eq. ( p.8[ ) into Eq. ( |1.9| ). Here, the eigenvalue A 
represents the multiplicative growth or decay of the total number of walkers that occurs at 
each time step. Since we are interested in distributions of random walkers for which the 
birth rate balances the death rate (so that the total number of walkers is constant in time), 
we must set A = 1 in Eq. (|2l^). We solve this eigenvalue equation for the case Z) = 1 in 



>^9n = < 



Sec. Ill and for the case of arbitrary D in Sec. [V 



III. ONE-DIMENSIONAL RANDOM WALKS WITH BIRTH AND DEATH 

In this section we consider the one-dimensional {D = 1) version of the discrete eigenvalue 
problem Eq. (|2]|). When D = 1, Eqs. ( iL^lLGl) and ([T^) both reduce to 
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Pnnt(n) 



i (n>2) 



1 (n = l] 



and 

in>2), 

Pinin) = I 

[o {n = l). 

For this case the steady-state distribution obtained by setting A = 1 in Eq.( p.9|) satisfies 

^zgn+i + \zgn-i (n > 3), 

\zg^ + agi in = 2), (3.1) 

. \zg2 {n = l). 

It is easy to solve the difference equation because it is a hnear constant-coefficient 
equation. Its general solution has the form 

= Ar"^-^ + Brl'^ {n>2), (3.2) 

where 

ri - -r± + 1 = (3.3) 

z 



and A and B are arbitrary constants. The solutions to the quadratic equation (|3.3| ) are 

1 



r± = - (l ± Vl - z^) . (3.4) 
Observe that 

r_r+ = 1. (3.5) 

Since the total number of random walkers is finite, the sum J2^=i 9n exists. From the 
existence of this sum and Eq. ( p.5| ) we may conclude that r±. are real; if r± were complex 
then, since they are complex conjugates, we would have |r-|-| = 1 and the sum would diverge. 
Furthermore, since r+ > 1, it follows that B = 0. 

If we substitute the solution ( ^.21 ) with B = into the special cases (n = 1 and n = 2) 
of Eq. (|3.1| ), we obtain a relationship between the birth rate a and the death rate z: 
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a = r I = - f 1 + Vl — , or z = — . (3.6) 

z ^ ^ a'' + 1 

This is the equation for the curved part of B2, the boundary curve between the upper region 
where Nt ^ 00 and the lower region where Nt ^ when a > 1. 

As a function of the birth rate a, the fraction F{a) of random walkers in region 1 for 
the steady-state distribution gn is given by 

= i = 

(Note that the overall multiplicative constant A drops out from this result and is unimpor- 
tant.) Equation ( |3.7| ) is only valid for a > 1; if a < 1 then no nontrivial steady-state solution 
exists; the limiting value of Ff as t — > 00 is 0. Indeed, it is shown in the Appendix that as 
t ^ 00 the fraction Ft vanishes like l/y/i along the line Bi and like 1/t everywhere to the 
left of Bi. 

We observe a second-order phase transition in F{a) = \imt-,oo Ft{a) as a function of the 
birth rate a; below the critical birth rate Oc = 1 this fraction vanishes and just above the 
critical point the fraction rises linearly with slope 1/2: 

F(a)~^(a-ae) (a ^ 1+). (3.8) 



Hence, at D = 1 the critical exponent u in Eq. ( |1.1CI| ) is 1 and the constant C{1 



2' 



IV. D-DIMENSIONAL RANDOM WALKS WITH BIRTH AND DEATH 

In this section we generalize the analysis of the previous section to arbitrary dimension 
D. When D ^ 1 the difference equation ( |2.9| ) is no longer a constant-coefficient difference 
equation and it cannot be solved in closed form. Hence, we use the method of generating 
functions to study steady-state (A = 1) solutions of this difference equation. 

We seek a solution to the D-dimensional generalization of Eq. ( ^.11 ) 
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9n — \ £)_,_3 
1 

D+1 



zg3 + agi 



zg2 



[n 



[n 



2) 
1) 



(4.1) 



which is obtained by substituting the uniform approximations to Pout(^) and Pmin) given 
in Eq. (^) into Eq. and setting A = 1. 

For a steady-state solution having a finite number of random walkers the sum 9n 
exists. We may therefore sum both sides of Eq. ( |4.1|) from n = 1 to oo and simplify the 
result: 



J29n = (a- z)gi + z^gr, 



(4.2) 



n=l 



n=l 



Assuming that the sum Y^^=i 9n is nonzero we may immediately conclude that 

E„=i gn a - z{a) 

Note that the result in Eq. ( [4.3|) is valid on the curved part of B2, where the sum exists and 
is nonzero; it is also valid on the straight-line portion of B2 when D > 2. On the curve B2 
we must treat 2; as a function of a. We emphasize this dependence by writing z{a) and by 
treating the fraction F as a function of a only. 

To obtain the dependence of ^ on a along B2 we use generating function methods. To 
begin, we simplify Eq. (|4. 1|) by setting 



gn = {2n + D - 3)/i„. 

Substituting Eq. ( ^^ into Eq. ( ^ ) we obtain 

' nzhn+i + {n + D — 3)zhn-i {n >3), 



(4.4) 



{2n + D- 3)K = { 2zh3 + {D - l)ah 



Zhq 



(n = 2), 
(n = l). 



(4.5) 



Next, we define a generating function: 



H{x) = Y.x'^h 



n+l- 



(4.6) 



n=0 
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Note that 



G{x) = x^gn+i = [2x— + D-lj H{x). 

Multiplying Eq. ( |4.5| ) by a;"~^ and summing both sides from n 
first-order inhomogeneous linear differential equation for H{x): 



(4.7) 



3 to cxD, we obtain a 



[zx^ -2x + z)H'{x) + {D- l){zx - l)H{x) = (z - a){D - l)xhi. 



(4.8) 



To solve Eq. ( [4.8|) we multiply both sides by the integrating factor {x'^z — 2x + z) 2 . The 
differential equation then simplifies to 
d 



dx 



zx — 2x + z) 2 H{x) = (z — a)gix{x z — 2x + z) 



D-d. 
2 



(4.9) 



The general solution to Eq. ( [4.9| ) is 



H{x) = [zx^ —2x + z 



l-D 
2 



C + (z - a)gi ds sis z -2s + z)~ 
Jo 



(4.10) 



where C is an arbitrary constant. 



To determine the constant C we observe that H{0) = hi = from which it follows 



that 



Thus, 



X 



H{x) = giix' -2- + 1 



1 



a 



. A--] I dss[s^-2- + l 
D-1 \ zj Jo V z 



D-3 

2 



(4.11) 



Finally, we use Eq. ( |4.7| ) to obtain the generating function G{x): 



X 



G{x) =gi< x'-2- + l {l-x') 



D+l 
2 



^ + («-^'ll-f)/>-(^=-2^ + l 



D-3 
2 



(1 - f ) 
x2-2f + 1 



(4.12) 



Assuming that G{x) exists for all < x < 1, we formally recover Eq. ( [4.3] ) when we set 
X = 1. 



14 



Recall the quantities r± defined in Eq. (|3.4|) and rewrite Eq. ( [4.12|) as 



G{x) = gi< [{r+ - x){r^ - x)] ^ (1 - a; ) 



X / dss{r+ — s) 2 (r_ — s) 



D-S 
2 



2a;2 ( 1 - f 



(4.13) 



(r+ — x)(r_ — x) J 

For z < 1 the generating function G{x) may be singulare at x = r_ < 1, in which case 
the representation of G{1) as a series will not exist. To preclude the possibility of such a 
singularity it is necessary and sufficient to impose the following eigenvalue condition: 

/„ ^ / Cb\ f^- , , , -D-3 , , D-3 

1 + (D - 1) f 1 - -j y dss{r+-s) — {r_-s)—=0. (4.14) 

This condition is clearly necessary. We can verify that it is sufficient by showing that 
G{r_ — e) exists in the limit as e ^ 0"''. To leading order in e the eigenvalue condition in 
Eq. ( [4.14|) becomes 

l + {D-l)h--j J dss{r+~ s) — {r^- s) — 

/ 0'\ , , -D-3 -D-1 ,1, , , 

~ -2 M - -j r_(r+ -r_) — e— (e ^ 0+). (4.15) 

Substituting this asymptotic result into Eq. ([4.13|) we see that the last term, which is of 
order e~^, exactly cancels. 

The eigenvalue condition in Eq. (|4.14| ) expresses the relation between a and z that we 
seek. We can rewrite this condition more compactly by rescaling the integration variable. 
Let s = r_u, so that 

{D - l)r^_ - 1^ duu{l - u)^{l-r^_u)^ = 1. (4.16) 

This integral converges only if D > 1. We can analytically continue to values < D < 1 
[recall that the region of validity of the uniform approximation in Eq. (|1 . 7| ) is D > 0] by 
recognizing that this expression contains the standard integral representation for a hyper- 
geometric function 

4r2 fa \ /3-D D + 3 ^\ 

- ' --l)2FA^—,2;^—;rl) =1. (4.17) 



D + l\z J \ 2 
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Using the transformation formulas (especially 15.3.26, 15.2.18, and 15.2.20 in Ref. 0]) for 
hypergeometric functions, this eigenvalue condition can be simplified: 



z\ n D + i 



(4.18) 



This implicit equation determines the curve B2 in the (a, z) plane. However, such a higher 
transcendental equation cannot be solved for 2; as a function of a in closed form. Thus, we 
perform an asymptotic analysis of this condition for z near 1. [As in the previous section, 
we find that for z ^ 1~ along B2 there is a transition at 2; = 1 from nontrivial steady-state 



solutions to trivial solutions of the walk equation (j^A 

To perform this analysis we let z = 1 — 77. We then use the following formula for the 
analytic continuation of a hypergeometric function: 

2Fi(a, b; c- = r(c)r(c-a-&) + b - c + 1-1 - () 

T[c — a)T{c — b) 

+ (1 - ^)c-a-. r(c)r(a + & -c) ^^ _ c - 6; c - a - 6 + 1; 1 - C). (4.19) 
V[a)V{b) 

Next, we substitute the first few terms in the Taylor series of a hypergeometric function 

Tic 



2Fi{a,b; c; C) 



r(a)r(6) 



g r(a + n)r(& + n) ^. 



n=0 



n\T(c + n) 



(4.20) 



to obtain 



an — 1 a — ar ri 
+ ^ + - 



D-i r 2 

1 +r? 

D-2\ 'A-D 



K 



2~D 



(4.21) 



where 



(4.22) 



which is valid near the critical point {ac,Zc = 1). Note that the value of Oc depends on D 
and must be determined by Eq. ([4.21|) . 



Our results are as follows. Leading-order asymptotic analysis for small rj gives the lo- 
cation of the critical point (a^^^c); the critical point lies at (1, 1) for < -D < 2 and at 
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{D — 1,1) for D >2. A next-order asymptotic analysis of Eq. ([4.21|) for the case < D < 2 
yields 



K \ 2-D 2 

2:(a) ~ 1 - ( ^— I (a-ar)~o [a ^ a; 



and 



F(a) 



K \ 2-D D 

(a — ar)^^ fa — > 



2-D 



Equation ( 4.24 ) reduces to Eq. (|3.8| ) when D = 1. 

For the case 2 < D < A a next-order asymptotic analysis of Eq. ( |4.21 ) gives 



z{a) r^l-[{D-2)K]T^ {a-a^)T^^ (a ^ 



and 



F{a) ^ [{D -2)K]~ {a- a^)— (a ^ a+). 



When D > 4 we find that 



D - A 

z{a) ^ 1 - (a-^flc 



and 



F{a) 



D-A 



-(a — ttc) (a ^ a^}"). 



1)(D- 2) 

The special case D = 3 can be solved exactly in closed form: 

fa -2)2 



zia) = 1 — 



o? + A 



{a > ac = 2) 



and 



F(a) 



(a - 2f 



fct ^ (^rmc 2) 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



(4.28) 



(4.29) 



(4.30) 



Indeed, the difference equation (|4.1| ) can be solved exactly and in closed form for all odd- 
integer D] the solution that vanishes as — * cxo is given by 

9n = r'lPD_i{n), 
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where Vk{n) is a polynomial in the variable n of degree k. When D is an odd integer the 
hypergeometric series in Eq. ( [4.17[ ) truncates for D > 5. Unfortunately, except for the cases 



D = 1 and D = 3 we do not obtain a simple form for the solution for z{a) and F{a). An 
implicit solution for z{a) when D = 5, for example, is given by 

{9a^ + 64)^3 - 56az2 - (Sa^ + 48)^ + 48a = 0. 

The special cases D = 0, D = 2, and D = 4 need to be treated separately. For D = 
the eigenvalue condition Eq. (|4.17| ) becomes very simple because we can use the identity 

2Fi(a,6;a;C) = (l-C)^'. 

Elementary algebra then yields 

z(a) = - (4.31) 
a 

for all a. Thus, the boundary B2 is a hyperbola for all a; the straight-line portion of B2 for 
a < 1 disappears. To understand this result observe that when D = Eq. ( [4.1| ) states that 
random walkers in region 2 cannot move outward. The appearance of this restriction is an 
artifact of the uniform approximation in Eq. ( 11.7] ). Thus, a steady-state solution has Qn = 
for n > 2 and consists of random walkers oscillating between region 1 and region 2. In this 
case, the fraction F{a) of walkers in region 1 is exactly 

F{a) = (4.32) 
1 + a 

For this degenerate case there is no critical point and no phase transition. We emphasize 
that the disappearance of a phase transition is an artifact; the uniform approximation in 
Eq. ( |1.7| ) is only valid when D > 0. 
For D = 2 we find that 



2 



z{a) ~ 1 — (constant) e (a ajT) (4.33) 

and 



For Z) = 4 we have 



z{a) ~ 1 



a — ttc 



and 



F{a) 



61og(^) 



a — ac 



a — > aj 



(4.35) 



(a ^ a+). 



(4.36) 



12 log 

The resuhs in Eqs. (|4.23| - |436|) confirm the formulas given in Eqs. (|1.10| - |rT2D . 

The limiting case D ^ oo is interesting because, like the case D = 1, we can find 
the exact equation for the curved portion of i?2- To treat this case we perform a large- 
asymptotic expansion of the integral in the the eigenvalue condition given in Eq. ( [4.161) . 
Using Laplace's method we obtain an asymptotic expansion of this condition as a formal 
series in powers of 1/Z^. We recover from this condition an expression for z as a. function of 



z{a) + — 

a D 



a 



(4.37) 



As one can see from Eq. ( [4.28] ) , in the limit D ^ oo the transition at a = Oc = -D — 1 is still 
second order. However, in this limit the discontinuity in the slope of F{a) disappears and 
F{a) for all a. 

Equations ( |4.25| ) and ( [4.27|) indicate that there is a change in the form of the transition 
at Z) = 4. When D < 4 the slope of the boundary curve B2 is continuous, and the critical 
exponent depends on D. However, when Z) > 4 an elbow appears in B2 at the critical value 
Oc = D — 1, and the critical exponent is independent of D. Specifically, when D > 4 the 
slope of Z?2 is for < a < Z^ — 1; just above a = Z^ — 1 the slope abruptly becomes 



D-4 



D{D-1)- 

We conclude this section by presenting a quick heuristic argument that reproduces the 
results in Eqs. ( [4.27| - [4.28D . For the case Z) > 4 we showed in Ref. that Ti{D), the expected 
time for a random walker who originates in region 1 to return to region 1, is given by 
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Tip) = 2^. (4.38) 

In a steady state all gi random walkers in region 1 leave this region and in a 2; = 1 model 
only the fraction Ili{D) ever return. The random walkers who return to region 1 do so 
in Ti{D) steps on average. These returning random walkers experience a death rate z for 
Ti{D) — 1 of these Ti{D) steps. Thus, the expected number of random walkers who actually 
return to region 1 is decreased by the factor z^^^^"*^^ . Hence, after Ti(D) steps we expect 
to find aIli{D) z^^^^^~^ gi random walkers in region 1. The condition that there be a steady 
state is therefore given by 

ani(D)z^^(^)-^ = 1. (4.39) 

Using the expressions for Ili{D) and Ti{D) in Eqs. ( ^.71 ) and ( ^.381 ) we obtain an approximate 
relation between z and a that is valid near the critical point; that is, where a = D — 1 + 6, 
z = 1 — e as (5, e — s> 0+. To first order in 6 and e this approximate relation is 

which is precisely the result in Eq. (|4.27| ). We obtain the result in Eq. ( [4.2^ ) by substituting 
Eq. ( [4. 271) into Eq. ( [4.3|) . Note that this argument is valid only for a > ac = D — 1. 

While the above argument is only valid in the neighborhood of ac, we can also use the 
above reasoning to derive the entire curve z{a) in the limit D 00. In this limit, Ti = 2. 
Hence, from Eq. ( [4.39| ) we have 

z = —, 
a 

the leading behavior in Eq. ( [4.37|) . 



(4.41) 
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APPENDIX A: LARGE-TIME ASYMPTOTIC BEHAVIOR 



In this Appendix we analyze the large-t behavior of the distribution Gn,t of random 
walkers for the uniform approximation of the probabilities Pout and Pin given in Eq. { \L.7\) . 



To this end we solve the set of equations in (L9) for the Kronecker delta initial condition 
Gn,o = Sn,i- As discussed earlier, the large-t behavior of a dissipative process is independent 
of the specific choice of initial condition. 

First, we derive a formal solution for Gn,t that is valid for general Pout and Pin- We define 



dn,t - 

and rewrite Eqs. (|1.9|) as 



az 



n-2 



out V 



n,t 



n>2), 
1), 



[n 



Qndn+l,t-l 



i-1 



n,t 



where we let 



Qid2,t-i 



Qn 



Next, we define the generating function 

oo 

en{y) = ^dn,ty^ 

t=0 



{n > 2), 
(n=l), 

(n > 2), 
(n = 1). 



and obtain from Eqs. (|A2| 

en{y) 



n>2) 
1) 



(n 



yQnen+i{y) +yen-i{y) 
l+yQie2{y) 

where we have applied the Kronecker delta initial condition. 
Let us define a continued fraction by the recursion relation 

1 



Sniy' 



n > 1 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



1 - y'^QnSn+iiy'^) 

It is easy to show that for n > 3 the recursion relation in Eq. (|A5|) is satisfied by 



(A6) 
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(A7) 



i=2 



[Since eniy) obeys a second-order difference equation, there is a linearly independent solution 
which can be determined using the technique of variation of parameters. This solution does 
not contribute; apparently, it fails to obey the appropriate boundary conditions at = oo.] 
We determine ei and the constant A by solving simultaneously the special cases n = 1 and 
n = 2 of Eqs. (|X5|): 



AyS,{y') = AQ,y'S2{y')Ss{y')+yeM, 
ei{y) = l + Ay'Q,S2{y'). 



(A8) 



Solving the above equations leads to a surprisingly compact expression for all en{y)- 

n 

en(l/)=l/"-' l[S^iy') {n>l). (A9) 



i=l 



Using a contour integral to project out the coefficients in the generating function we obtain 



G 



n,t 



z, 



n-l 



out V 



Li=l 



dy 
c 2my 



n-t-l 



(AlO) 



where empty products are defined to be unity. The contour C encircles the pole at the origin 
in the complex-?/ plane but excludes all other singularities of the integrand. 

This rather strange expression (a contour integral over a product of continued fractions!) 
is of little use, even in an asymptotic analysis for large values of t. Only for particular choices 
for Pout and Pin is progress possible. A significant advantage of the uniform approximation 
in Eq. ( |1 . 7|) [compared with probabilities given in Eqs. ( L4| - T^ )] is that they simplify the 
expression for Gn,t in Eq. ( AlO ) for all D > 0. [The probabilities in Eqs. (|1.4| - |1.6|) ] lead to 
a tractable result only for D = 0, 1, and 2.] 



We simplify the expression for Gn,t in Eq. (|A10|) by recalling the continued fraction 



representation for a hypergeometric function |[TT 

2Pi(a,6 + l;c+l;C) 



2Pi(a,6; c; C) 



1/(1 + /iC/(l + /2C/(l + /3C/(l + ...)))), 



(All) 
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where 



{i + b){i + c-a) 



2i 



and /si+i 



{i + a){i + c — a) 



{2i + c){2i + c-l) {2i + c){2i + c + 1)' 

Substituting the uniform approximation for Pout and Pin in Eq. ( [1.7| ) into Eq. 

n{n + D-2) 



Qn ^ 

which can be rewritten as 

Qn+2i = z'^ 



{2n + D -3){2n + D -1) 



{n > 2), 



(A12) 
gives 
(A13) 



n+D-2 



)(^ + f) 



Q 



n+2i+l 



{2i + n + ^){2i + n + ^^ 



in > 2A > 11 



n nTTT (n>2,i>0). 

{2t + n + ^){2i + n + ^) ^ - ' - ^ 



(AM) 



The continued fractions in Eq. ( |Aq ) can thus be identified as 

P (n n+1 I D-1 . 2 2 



g-3 . ^2„2 

2 ' ^ y 



(n > 2). 



2Fi(f,V;n + 

Hypergeometric functions are symmetric in their first two arguments. Therefore, 

2Fi{l^;n + ^;zV^ 



(A15) 



i=2 



2-^1 (^2' ^ y ) 

Substituting this last resuh into Eq. (|A10|) we finally obtain 



{n > 2). 



(A16) 



Gn.t — ^ 



2"-2r(D)r (n + ^) 2my ^ 1 + (f - l) 2^^! (i 1 



2 ' 2 1 " ~r 2 ' y 



(A17) 



Recalling Eq. we obtain an expression for the total number of walkers at time t by 
summing Eqs. ( p..9|) over all positive integers n: 



t-i 



T=0 



(A18) 



Next we insert Gi^t from Eq. ( [A17| ) and sum over r: 



AT, = 2* 



1+ 1-- 



2^1 I 2' 1' 



(A19) 



23 



where we have ehminated terms in the integrand that are regular at the origin in the complex- 
y plane. From Gi j in Eq. ( |A17[ ) and Nt in Eq. ( |A19| ) we obtain the large-t behavior of the 
fraction Ft in Eq. ( p.2| ). Note the similarity of the denominator in both integrals with 
the eigenvalue condition in ( [4. 181 ) for steady-state solutions. The asymptotic behavior of 
the integrals for large t is dominated by the poles of the integrands, and the steady-state 
solution is merely the special case where the asymptotic behavior is independent of t in 
leading order. 

To extract the large-t behavior of Nt and Ft, we conduct a saddle-point analysis of the 
contour integrals for Gi^t and Nt in Eqs. ([A17|) and ( |A19|) . Both expressions, aside from the 
prefactor z*, only depend on the ratio a/z. Saddle-point analysis requires that we consider 
three distinct cases, values of a/z such that (i) (a, z) lies to the left of the line (ii) (a, z) 
is on and (iii) (a, z) lies to the right of Bi (see Fig. [l|). For all cases the integrands for 
both Gi^t and Nt have a pole aX y = on the real positive axis. For case (i) t/p > 1, for 
case (ii) ?/p = 1, and for case (iii) y^ <1. 

For case (i) we find to leading order that as t ^ oo 
- {2-DY 



Git ~ < 



2KT 



zH"^-' 



(D-l)z 



2^-\D-2fV 



D-l 



{D-l)z 



TT 



(0 < L> < 2), 
{D > 2), 



(A20) 



and 



{2-D) , D . 



1 - f ) irr (f ) 



iD-l)z 



1 - 



-{D - 2)z' 



{0< D < 2), 
{D > 2), 



(A21) 



{D-l)z 

where K is given in Eq. ( 4.22| ). Hence, 



r 1 



Ft ~ < 



D 



2T-2(/)_2)r 



D-l 



{D-l)z 



D 



TT 



(0 < D < 2) 
{D > 2). 



(A22) 
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For case (ii) we find to leading order that as i — > oo 

K 



Git ~ ' 



2r(2-f 



D-l , D ^ 



2KT (f 

D-A 
I 2{D-2)' 



and 



D-l 



Nt ~ < 



KT 



D 



D-A 



Hence, 



Ft ~ <^ 



-zH 



D-l D 



ifr(f) 

1 

2t 

1 



(0 < D < 2), 
(2 < L> < 4), 
{D > 4), 

(0 < L> < 2), 
{2<D < 4), 

(D > 4). 

(0 < D< 2), 

(2 <D <4), 
{D > 4). 



(A23) 



(A24) 



(A25) 



I {D - 2)t 

The analysis of case (iii) is somewhat more comphcated. The saddle point in cases 
(i) and (ii) is very near y = 1 for large but in case (iii) the integrands have poles at 
Q < y — yp{a/z) < 1 and the saddle point is now located near yp. An asymptotic analysis 
of this case for large t is possible only if we consider a small neighborhood to the right of 
the hne Bi. Approaching Bi we find that yp — > 1". We use the ansatz ^ — ^ + e and 
yp{a/z) = 1 — 5(e), where e -C 1 and 5 <C 1, but where t5 ^ 1. We find that as i ^ oo 

K 



2 

£2-13 



2 

2-D 



2-D\ 

[K{D - 2)]"^ 

D-A 
{'2{D-l){D-2) 



{0<D< 2), 
(2<D< A), 
(D > A), 



(A26) 



leading to 
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Gi,t 



D 
£ 2-D 



K 



2-D 
4^ 2{D-1 



2 

2-D 



D-2 



2/p(aA) 



D-A 
2{D-2) 



1 1 



ypia/z)_ 



{0<D < 2), 
{2<D < 4), 
(D>4) 



and 



D 



2{D - 1) 



?/p(a/x 



e-^ (Z^-l)(D-2) 



(0 < D < 2), 
(2 < D < 4), 
(D > 4). 



Hence, 



D 
g 2-13 



K 



2 

2-D 



2-D. 
\K{D - 2)\ 



D-2 
D-A 



{0<D < 2), 

{2<D < 4), 
{D > 4). 



(A27) 



(A28) 



(A29) 



2{D-l){D-2y 

From the previous formula we can recover the asymptotic results in Eqs. (|4.24| ), ( [4.26[ ), 
and ([4.28|) , which are valid on the line B2 as a ^ a+ and z ^ z~ . The particular path B2 is 
distinguished merely by the fact that the total number of walkers Nf approaches a nonzero 
constant as t ^ 00. Thus, the portion of B2 to the right of Bi is obtained for z = yp{a/z) 
in Eqs. ( |A28| ). Again, we let 2; = 1 — 77 for r/ 0"*" and find that rj ~ S{e). Then, using 
Zc = 1, we find that e ~ (a — ac) + acTj. From Eq. ( |A26| ) for D < 4, we have e ^ rj, and we 
merely need to identify e = a — Oc in Eq. ( [A29| ) to recover our earlier results. For D > 4 we 
recall that e = 0{r]) to recover Eq. ( |4.28| ). 
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FIGURES 

FIG. 1. Generic phase diagram for the (a, z) plane. (The diagram was actually generated using 
data from D = ^ random walks.) Shown on the diagram are the boundary curves Bi and B2. To 
the left of Bi and on Bi the fraction of random walkers in region 1, Ft, approaches as t ^ 00; 
to the right of Bi this fraction approaches a finite positive number as t ^ 00. Above B2 the 
total number of random walkers, Nf diverges as t ^ 00; below B2 the total number of walkers 
approaches as t ^ 00. On B2 the distribution of random walkers approaches a steady state as 
t — > 00. The critical point (oc, Zc) lies at the intersection of Bi and i?2- 

FIG. 2. Phase diagram in the (a, z) plane for the case D = 1. For this dimension the proba- 
bilities Pout('^) and Pin{n) for the hyper spherical surface area case given in Eqs. (1.4-1^) and the 



uniform approximation case given in Eqs. (L7) are the same. On this diagram the slope of Bi is 



unity, ttc = 1, and the slope of B2 is continuous. 

FIG. 3. Phase diagram in the (a, z) plane for the case D = 2. For this dimension the probabili- 
ties i-'out('^) and Pin{n) for the hyperspherical surface area case and the uniform approximation case 
are the same. On this diagram the slope of Bi is unity, Oc = 1, and the slope of B2 is continuous. 
Note that the slope of B2 approaches exponentially fast as a — > Oc from above. 

FIG. 4. Phase diagram in the {a,z) plane for the case D = 3 using the probabilities Pouti'n) 
and Pin{n) for the uniform approximation case in Eqs. (|1.7D . The slope of Bi is 1/2, Oc = 2, and 
the slope of B2 is continuous. 

FIG. 5. Phase diagram in the (a, z) plane for the case D = 2> using the probabilities Pout("') 
and Piri{n) for the hyperspherical surface area case given in Eqs. ( |1. 4^1.61 ). The slope of Bi is 
1 — 1/C(2) = 1 — 6/7r^, Oc = 2.551 . . ., and the slope of B2 is continuous. 

FIG. 6. Phase diagram in the (a, z) plane for the case D = A using the probabilities -Pout('^) and 
Phi{n) for the uniform approximation case given in Eqs. ( [1.7] ). The slope of Bi is 1/3 and Cc = 3. 
The slope of B2 is continuous; it vanishes logarithmically as a ^ ac from above [see Eq. ( [4.35 )]. 
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FIG. 7. Phase diagram in the (a, z) plane for the case I? = 4 using the probabilities Pout('T') 
and Pinin) for the hyper spherical surface area case given in Eqs. ( |1.4Hl.6| ). The slope of Bi is 
1 — 1/C(3) and ac = 5.949. . .. Universality arguments lead us to believe that the slope of B2 is 
continuous and vanishes logarithmically as a — > Cc from above as in Fig. |^. 

FIG. 8. Phase diagram in the {a,z) plane for the case D = 5 using the probabilities -Pout("') 
and Pm{n) for the uniform approximation case given in Eqs. ( |1.7D . The slope of i?i is 1/4 and 
Oc = 4. The slope of B2 is not continuous; there is an elbow at a = Oc- 

FIG. 9. Phase diagram in the {a,z) plane for the case D = 5 using the probabilities Pout{n') 
and Pinin) for the hyper spherical surface area case given in Eqs. (1.4-L^). The slope of Bi is 



1 — 1/C(4) = 1 — 90/7r^ and ac = 13.147. . .. The slope of B2 is not continuous; there is an elbow 
at a = Qc- 
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